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ABSTRACT 

We present results of a 2D3V kinetic Vlasov simulation of the Weibel instability. The kinetic Vlasov 
simulation allows us to investigate the velocity distribution of dilute plasmas, in which the effect of 
collisions between particles is negligible, and has the advantage that the accuracy of the calculated 
velocity distribution does not depend on the density of plasmas at each point in the physical space. 
We succeed in reproducing some features of the Weibel instability shown by other simulations, for 
example, the exponentially growing phase, the saturation of the magnetic field strength, the formation 
of filamentary structure, and the coalescence of the filaments. Especially, we concentrate on the 
behavior of the filaments after the saturation of the magnetic field strength and find that there 
is a kind of quasi-equilibrium states before the coalescence occurs. Furthermore, it is found that 
an analytical solution for stationary states of the 2D3V Vlasov-Maxwell system can reproduce some 
dominant features of the quasi-equilibrium, e.g, the configuration of the magnetic field and the velocity 
distribution at each point. The analytical expression could give a plausible model for the transition 
layer of a collisionless shock where a strong magnetic field generated by the Weibel instability provides 
an effective dissipation process instead of collisions between particles. 

Subject headings: plasmas — instabilities — shock waves — magnetic fields 



1. INTRODUCTION 

It has been a long time since the existence of collision- 
less shocks was confirmed. In this context, the collision- 
less shock means discontinuities of physical quantities oc- 
curring in sufficiently rarefied plasmas, in which the effect 
of collisions between particles composing the plasma is 
negligible. So far, it is believed that collisionless shocks 
play important roles in a wide class of astrophysical phe- 
nomena, e.g., supernova remnants(SNRs), the prompt 
emission and afterglow of gamma-ray bursts(GRB), jets 
from active galactic nuclei. 

From the discovery of t he bow shock in front of the 
terrestrial magnetosphere (|Ness et all 1 1 9641 ). the proper- 
ties of collisionless shocks, e.g., the mechanism of the 
formation, the shock conditions, and the configuration 
of the magnetic field, hav e been inv e stigat ed by various 
approaches. For example. iPaul et al (|1965T ) observed the 
formation of a collisionless shock in laboratory experi- 
ments. Observationally, several satellites have revealed 
the nature of collisionless shocks. Especially, GEOTAIL 
measured the variation of the magnetic field strength in 
the bow shock in front of the terrestrial magnetosphere 
from the upstream to the downstream. GEOTAIL also 
acquired the energy spectra of ions in the solar wind, 
which confirmed the existence of non-thermal particles. 

Also in other high-energy astrophysical phenomena 
such as SNRs and GRBs, the existence of non-thermal 
particles is confirmed observationally. In generating such 
particles, a strong magnetic field amplified around the 
collisionless shock is considered to be a key ingredient. 
This prospect is supported by observations of several 
SNRs tiVink fc Laming! 12001 iBamba eFaXl 120031 ). Ad- 
ditionally, the recent observation of the time variation 

1 Research Center for the Early Universe, School of Science, Uni- 
versity of Tokyo, Bunkyo-ku, Tokyo 113-0033, Japan. 

2 Department of Astronomy, School of Science, University of 
Tokyo, Bunkyo-ku, Tokyo 113-0033, Japan. 



of X-r ay sources in RX J1713. 7-3946 bv lUchivama et all 
(|2007h strongly indicates the existence of a strong mag- 
netic field in the shock front. Also in GRB afterglows, 
the existence of a long-lived, near equipartition magnetic 
field is indisp ensable for the observed syn chrotron emis- 
sions(see, e.g, lRrar] [l999: Mcs zlr^[2006l) . 

Theoretically, a great deal of investigation has been 
done in both analytical and numerical framework. 
iMedvedev and Loeb (1999) s uggest ed that the Weibel 
instability(Weibel 1959c iFriedl 119591 ). which is a kind of 
plasma instability caused by the anisotropic velocity dis- 
tributions of collisionless plasmas, could dissipate the ki- 
netic energy of the plasma and generate a strong mag- 
netic field in the transition layer of a collisionless shock. 
First principle simulations of the formation of a colli- 
sionless shock (jKatol 120071 : ISpitkovskvl l2008al|bl ) by us- 
ing recent powerful computers have su pported the sug- 
gestion. Especially, ISpitkovskvl (|2008bD performed long- 
term 2D particle-in-cell(PIC) simulations of a relativistic 
collisionless shock in an unmagnetized pair plasma and 
found a sign of the generation of non-thermal p articles 
via t he first-order Fermi acceleration mechanism (jFermil 
1949). They conclude that the scattering source of parti- 
cles, which is needed for the Fermi acceleration to occur, 
is a turbulent magnetic field generated by the Weibel in- 
stability in the transition layer of the shock. This is why 
the Weibel instability is a compelling candidate for the 
mechanism of the formation of the shock front associated 
with a strong magnetic field and particle acceleration. 

However, there remain some problems regarding the 
Weibel instability itself. Results o f recent 3D sim- 
ulations of the Weib el instability (jSilva et alj 120031 : 
iFrederiksen et al1l2004D show that filamentary structure 
(sometimes called "current filament") forms as a mag- 
netic field is amplified. After saturation of the magnetic 
field strength, the filaments coalesce each other to form 
larger one. Although each filament is small (its initial 
radius is roughly equal to a few electron skin depth), 
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the detailed properties of this structure, e.g., the cor- 
relation length of the magnetic field, the particle en- 
ergy distribution, and the time evolution, may affect the 
macroscopic behaviors of the transition l ayer o f a colli- 
sionless shock. For example, iMedvedevl 1)20001 ) pointed 
out a possibility that ultrarelativistic electrons acceler- 
ated by highly nonuniform small scale magnetic fields 
emit radiations different from a synchrotron radiation 
emitted by electrons in uniform magnetic fields. Subse- 
quently some investigations into the filamentary struc- 
ture have been perfor med. From the conservation of en- 
ergy, iGruzinovl (|2001h pr edicted that t he ma gnetic field 
energy decreases as i -1 . IChang et al.l (|2008l ) considered 
a variant of Landau damping (a process dissipating the 
energy of electromagnetic fields into the kinetic energy 
of particles) and reached a similar result for the evolu- 
tion of the magnetic fields. They assumed straight or- 
bits of particles because of weak magnetic fields, though 
their PIC simulations reveal that there exist some re- 
gions whe re strong magne t ic fiel ds bend the orbits of 
particles. M edvedev et al.l (|2005l ) investigated the coa- 
lescence of the filaments by using a two-dimensional toy 
model. They pointed out that the correlation length of 
the magnetic field g rows exponentially at fir st and then 
linearly with time. lAchterberg et alT (|2007f ) reexamine 
the toy model including the effect of screening currents of 
the background plasma and predict that the coalescence 
slows down when the separation of the filaments becomes 
larger than the skin depth of the background plasma. 
But their analysis does not predict whether there exists a 
critical scale length above which the coalescence does not 
proceed. Furthermore, their simple model ignores micro 
processes such as the dissipation of the magnetic energy 
into the kinetic energy by the re connection of the mag- 
netic field. On the other hand, iMilosavliev ic fc Nakerl 
(2006) treated the filamentary structure in the frame- 
work of the magnetohydrodynamics(MHD). They argue 
that a pressure-driven MHD instability destroys the fil- 
ament and the strong magnetic field generated in the 
transition layer is short-lived. That is, the consensus is 
yet to be reached regarding the long-term evolution of 
the Weibel filaments. Especially, it is unclear whether 
the toy model and the MHD model are appropriate to 
model the filaments. 

To settle such unresolved problems, we use a kinetic 
Vlasov code. This approach enables us to investigate 
the structure of the filaments and their coalescence and 
has the advantage that the accuracy of the calculated 
velocity distribution does not depend on the density of 
plasmas at each point in the physical space, while it does 
in PIC simulations. However, there are some computa- 
tional constrains on Vlasov codes. The most serious one 
is its small dynamic range in both physical and velocity 
spaces. Due to the defect, we cannot treat the formation 
of a shock from two plasma flows or relativistic motions 
of particles. Thus, in this work, we concentrate on the 
coalescence of filaments forming as a result of the non- 
relativistic Weibel instability. The saturation of the rel- 
ativistic Weibel instability is realized by the same mech- 
anism as that of the non-relativistic Weibel instability 
dKato fc Takabel[2008h . Therefore, if we can construct 
a reliable model of the filaments by analyzing the non- 
linear behavior of the non-relativistic Weibel instability 
in detail, it may provide a general insight into some es- 



sential processes operating in the transition layer of a 
collisionless shock and predict their long-term evolution 
beyond the reach of the simulation. 

This paper is organized as follows. In the next section, 
we describe the strategy for the numerical simulation. 
The results are shown in §3. In §4 , we construct an 
analytic model for the filaments and address remaining 
problems. We conclude this paper in §5. 

2. FORMULATION 

In this section, we describe the governing equations, 
some assumptions, and an initial setup for the numerical 
simulation. 

2.1. Equations 

The equations describing the behavior of rare fied plas- 
mas are the Vlasov-Maxwell system (see, e.g.. ISturrockl 
Il994f h The Vlasov equation governs the time evolution of 
the distribution function fj(t, x, y, z, v x , v y , v z ) of species 
j(i for ions and e for electrons). Here the cartesian coor- 
dinates are (x, y, z) and the corresponding coordinates in 
the velocity space are ( v x , v y , v z ) . The Maxwell equations 
govern the time evolution of the electromagnetic fields E 
and B. We impose two assumptions. One is that the 
plasma is homogeneous in the z direction, which makes 
the Vlasov equation to take the following form; 



dt 



' dx 



dv 



2f = 0, (1) 



dy in,-) 

where qj and rrij are the charge and the mass of species 
j, and c is the speed of light. E and B are expressed by 
introducing the scalar and the vector potentials, cj>(t, x, y) 
and A(i, x, y) as 



E 



-V4> =r, 

c at 



B = V x A. 



(2) 



The time evolution of the potentials is described by the 
wave equations; 

l01 " - 2 = -47rp, l5A-V 2 A = -4.j, (3) 



dt 2 



dt 2 



where p and j are the electric density and the electric 
current density, which satisfy the Lorenz condition; 



c at 



(4) 



The other assumption is that the plasma consists of ions 
and electrons with the same charge but with the opposite 
sign (qi = —q e = e). So the source terms in Equation (J3J) 
are expressed in terms of fj as 
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v(/< - f e )dv x dv y dv z . (6) 



2.2. Normalization 

We define some values characterizing the physical 
quantities; l/u> c as the time scale, c/w e as the length 
scale, c as the velocity, and m c cw c /e as the electromag- 
netic field. Here Lu e is the electron plasma frequency. 
Normalized by these values, the governing equations de- 
scribed in the previous subsection can be expressed in 
a non-dimensional form. In the following, we treat thus 
normalized physical variables. 
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2.3. Method of numerical integration 

To simulate the Weibel instability, we use the numeri- 
cal scheme for t he integration of the Vl asov-Maxwell sys- 
tem proposed by |MMiiniyI5i3 (EQQI). We review their 
procedure that evolves the distribution functions and the 
electromagnetic potentials by a time interval At in the 
following steps. 

In integrating the Vlasov equation (JXJ) , we treat the 
electromagnetic field as static. One can see the Vlasov 
equation |T]) consists of five one-dimensional advection 
equations in the form of 

where £(= x,y,v x ,v y ,v z ) is an independent variable, 
is a constant and g is a function of (t, (). This equation 
is solved by second-order van Leer's scheme. Here we 
define an operator T that evolves the function g by a 
time interval At according to Equation |7| ; 

g(t + At, t-a c At)=T[C At]g(t,Q. (8) 

Using this operator. iMangenev et all ()2002l ) proposed the 
following scheme to integrate the Vlasov equation (p}; 

fj(t + At, x, v) = S[x, y, At/2]S[v x ,v v , At/2] 

xT[z,At]S[v x ,v y ,At/2] 

xS[x,y,At/2]fj(t,x,v). (9) 

where another operator S has been defined as 

S[x, y, At] = T[x, At/2]T[y, At]T[x, At/2}. (10) 

If the distribution function at t + At is given, the 
source terms in Equations ([3]) are calculated by Equa- 
tion ([6]) . The wave equations ([3|) governing the time evo- 
lution of the electromagnetic fields are solved by using 
the Fourier transformation, which imposes the periodic 
boundary conditions in the x and y directions implicitly; 

/L x /2 rLy/2 
/ exp[-i(k x x + k y y)} 
-L x /2 J-Ly/2 

xh(t,x,y)dxdy, (11) 

where L x (L y ) is the length of the computational domain 
in the x(y) direction and h is a function of (t, x, y). The 
inverse Fourier transformation is defined by 

I rL x /2 rL y /2 

h(t,x,y)= / / exp[i(fc x x + k y y)] 

(2iryL x L y J-L v /2 

xh(t,k x ,k y )dxdy. (12) 

These transformations are computed by th e standard 
Fast- Fourier- Transformation scheme (see, e.g. lPress et all 
l2007f ) numerically. Substitution of these relations for 
(j),A,p,j into Equation leads to the following four 
second-order ordinary differential equations; 



rft 2 



-(k 2 x + k 2 v )^-4iTp, 



d 2 A 
~~aW 



+ A -47rj, 
(13) 

which are solved by the Runge-Kutta method of order 
4. To ensure that the obtained </>(t) and A(t) satisfy the 

Lorenz condition (H]), we use the values of d<j>/dt calcu- 
lated from (0} with A(t) instead of using those obtained 
by numerically integrating equations (|13p . 

In that way, we evolve the velocity distributions and 
the electromagnetic field one after the other. 



2.4. Simulation setups 

As the initial condition, we consider a counter- 
streaming plasma that is homogeneous in space; 



fjo = 



27r3/2 v 3 



exp 



- exp 



v x + vl + (v z -Vjf 



(14) 



where rij , Vj , and Vj are the number density, the thermal 
velocity, and the bulk velocity of species j. We assume 
that they are constants and n\ = n c — no for the charge 
neutrality. In our simulation, the values of these parame- 
ters are as follows; n^ = 1, v e = 0.05, Vi — 0.05y/m c /mi, 
and V c — Vi — 0.2. The mass ratio is assumed to be 
mi/m e = 1 and 16. The initial configuration of the elec- 
tromagnetic field is 

<f> = A x = A y = 0, A z = — e since sin y, (15) 
which is equivalent to 
E x = E y = E z = 0, 

B x = e sin a; cosy, B y — — e cos a; sin y, B z = 0(,16) 

where we have introduced a small parameter e(= 10 -5 ). 
In other words, we treat the above magnetic field as a 
perturbation to the initial distribution of particles with 
no electromagnetic field. Next, we address the boundary 
conditions. The simulation domain consists of the spatial 
intervals given by x, y G [— w, it] and the velocity intervals 
given by v x ,v y € [-0.4,0.4] and v z e [-0.6,0.6]. The 
periodic boundary condition is imposed in the x and y 
direction, while in the velocity space, the distribution 
function assumed to vanish for |t>a.|, \v y \ > 0.4 and \v z \ > 
0.6. The number of zones in the physical space is 32 x 32 
and that in the velocity space is 40 x 40 x 60. 

3. RESULTS 

In this section, we show the results of the integration 
of Equations (HJ)-® under the initial condition given in 
the previous section. 

3.1. The time evolution of energies 

Figure [1] shows the time evolution (0 < t < 200 for the 
case TOi/m G = 16) of the electron kinetic energies in the 
x and z directions, the electric energy, and the magnetic 
energy. The kinetic energy of species j in each direction 
is defined by 



7T />7T pOO />00 pOO 



K i*—2 III I I ' 

^ ./— 7T J —71 J —OO J —OO J— CO 



* J — 7T J — 7T J — OO J — OO J — OO 



%=af r r r r~< 

* J TV J 7T J — OO J — OO J OO 



xdxdydv x dv y dv z , (17) 
xdxdydv x dv y dv z , (18) 



xdxdydv x dv y dv z , (19) 
while the electric and the magnetic energies are defined 

by 

Ec = \j J ( E x+E 2 y + E 2 z )dxdy, (20) 
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Fig. 1. — The time evolution of the electron kinetic energies, K ex 
(solid), K ez (dashed), the electric energy E c (dash-dotted), and the 
magnetic energy E m (dotted), for < t < 200 and m,i/m e = 16. 
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Fig. 2. — The time evolution of the ion kinetic energies, Ki x 
(solid), Ki z (dashed), the electric energy E c (dash-dotted), and the 
magnetic energy E m (dotted), for < t < 600 and m,i/m e = 16. 
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{Bl + B 2 y + B 2 z )dxdy. 



(21) 



The values of Kj x and Kj y evolve in the same way, 
because of the symmetry in the v x -v y plane. The time 
evolution is divided into two phases, the linear phase 
and the saturation phase. In the linear phase, the 
electric and magnetic energies increase exponentially. 
The growth rate is c onsistent with a linearized analysis 
(|Califano et al.lll998l ). After the linear phase, non-linear 
effects become significant and the magnetic energy be- 
comes comparable to the total kinetic energy of electrons, 
while the electric energy takes relatively small values. In 
the transition between the two phases, K ex increases and 
K ez decreases. This is one of the remarkable features of 
the Weibel instability, the anisotropy of the motions is 
mediated by interactions between particles and electro- 
magnetic fields. 

From the time evolution of the ion kinetic energy 
in each direction and the electric and magnetic ener- 
gies shown in Figure [21 it is found that its large mass 
(mi = 16m e ) keeps the value of Ki Z almost constant by 
inhibiting the efficient exchange between the kinetic en- 
ergy and the electromagnetic energy. 

Figure [3] shows the time evolution (0 < t < 600) of 
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Fig. 3.— Same as Figured] but for < t < 600. 



K ex , K ez , E c , and E m . In the saturation phase, the 
electric energy gradually decreases indicating that the 
difference between the density distributions of ions and 
electrons diminishes. However, it does not mean that 
the contribution of the electric force to the force balance 
is negligible, because the balance of the Lorentz force 
acting on a charged particle is governed by the following 
relation; 

E + (Vj) x B = 0, (22) 
where (vj) is the bulk velocity of species j defined by 

oo />oo />oo 

/ / w] 1 dv x dv y dv z . (23) 

— oo J ~ oo J— oo 

In this simulation, a small value of the bulk velocity of 
ions of the order of |(vi)| ~ 0.1 enables the electric field 
to contribute to the achievement of the force balance. 

3.2. The saturated value of the magnetic field strength 

Several authors ( Califa ncTet al.l 119981 : 

iMedvedev and Loerj|l999t lKatdl2005f T have investigated 
the saturation mechanism of the Weibel instability. 
Their scenario is as follows. If the orbits of particles are 
straight, the particles can form current filaments, which 
generate magnetic fields around themselves. Thus the 
magnetic field strength grows consistently as long as 
the orbits of particles can be regarded as straight one. 
On the other hand, once the magnetic field strength 
becomes strong enough to deflect the orbits of particles, 
the particles cannot form current filaments. Hence, 
it is expected that the magnetic field strength will 
be saturated when the Larmor radius r^ of a particle 
becomes comparabl e to the chara c teristi c scale length 
I of the plasma. ICalifano et al.l (119981) adopted the 
electro n skin depth as I, while IMedvedev and Loebl 
(1999) adopted the reciprocal of the wave number of 
th e most unsta ble mode of the Weibel instability as /. 
In iKatol (|2005l ). results of a series of PIC simulations 
imply that the radius of a filament is compatible to the 
characteristic scale length. In our simulation, the wave 
number of the growing mode is fixed ab initio from 
the form of the perturbation (|15[) and identical to the 
electron skin depth c/w c , which leads to the formation 
of filaments with radius R = ir/2. Then, we obtain the 
following condition for the saturation of the magnetic 
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Fig. 4. — Snapshots of the time evolution of the Weibel instability. Each panel represents a color-coded density distribution of electrons (or 
ions) and the magnetic fields by arrows in the x-y plane. 
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We can estimate the saturated value of the magnetic 
field strength i? sa t = 2v c /ir = 0.13 by substituting 
V e = 0.2. The corresponding magnetic field energy 
is E m = Bl &t L x L v /2 = 0.32, while E m ~ 0.4 for 
nii = 16m and E m ~ 0.2 for m\ = m e from the re- 
sults of our simulation. Therefore, this rough estimation 
of the saturated value of the magnetic field strength is in 
good agreement with the results of the simulation. 

3.3. The density distributions 

Figure 2] shows snapshots of the electron (or ion) den- 
sity distribution and the magnetic field configuration at 
t = 90,200,300,340,380,400,420,460,500. These trace 
the time evolution of the filamentary structure resulting 
from the development of the Weibel instability. In the 
linear phase, the filamentary structure develops accord- 
ing to the initial configuration of the magnetic field. In 
the beginning of the saturation phase (90 < t < 300) , we 
can see two filaments carrying electric currents parallel 



to the z-axis around (x,y) = (— tt/2, — tt/2), (tt/2, tt/2) 
and the two other filaments carrying anti-parallel cur- 
rents around (x,y) — (n/2,—n/2), (— n/2, n/2). In this 
period, a quasi-equilibrium configuration seems to be 
achieved, since the distribution function hardly changes 
seems to be achieved. After a while, the configuration 
changes in the following way. Adjacent two filaments 
carrying currents to the same direction begin to ap- 
proach each other (300 < t < 340) and coalesce into 
a large filament (340 < t < 350). As the snapshots at 
t = 400, 420, 460, 500 exhibit, the large filament oscillates 
after the coalescence. The amplitude of the oscillation 
gradually decreases to achieve a new equilibrium config- 
uration, which is self-similar to the previous one. 

3.4. Thermalization of electrons 

As we see in the previous subsections, the Weibel in- 
stability amplifies a transverse electromagnetic field by 
rapidly dissipating the bulk kinetic energy of particles. 
However, all the dissipated energy is not converted into 
the electromagnetic energy. It is known that the electro- 
magnetic field thus amplified efficiently thcrmalizc elec- 
trons. We can see this effect in Figure where each line 
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Fig. 5. — The x, y, v y , ^-integrated electron velocity distribution 
at t = 10, 70, 100, and 300. 
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Fig. 7. — The ion velocity distribution at X-point (solid) and 
O-point (dashed) in arbitrary units for mi/m c = 16 and t = 280. 
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O-point (dashed) in arbitrary units for m\/m c = 16 and t = 280. 



represents the x, y, v y , w z -integrated velocity distribution 
of electrons at t = 10, 70, 100, and 300. Because the ve- 
locity distributions in Figure [5] are well represented by 
the Maxwell-Boltzmann distribution with zero bulk ve- 
locity, the kinetic energies K ex in Figures [1] and [3] are 
dominated by that of the thermal velocity. 
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Fig. 8. — The electron(or ion) velocity distribution in the quasi- 
equilibrium at X-point (solid) and O-point (dashed) in arbitrary 
units for m;/r?i c = 1. From the mass symmetry, the ion veloc- 
ity distribution at any point become identical with that of electron 
at the same point by transforming v% as v z — > — v z 



m; = m c . Because of the mass symmetry, the ion veloc- 
ity distribution becomes identical with that of electrons 
by transforming v z as v z — > — v z . The shape of the dis- 
tribution in Figure [8] is similar to that in Figure [6l 



3.5. The velocity distributions in the quasi- equilibrium 

In the first quasi-equilibrium (200 < t < 300), we 
focus on the velocity distribution at two points. They 
are the point at the center of a filament located at 
(x,y) — (7r/2,7r/2) (referred to as O-point), and the 
point surrounded by four filaments at (x, y) = (0, 0) (re- 
ferred to as X-point). In Figure solid and dashed lines 
represent the v x , i^-integrated electron velocity distribu- 
tions in the period of the quasi-equilibrium at the O- 
and X-points, respectively. At the O-point, the veloc- 
ity distribution is close to one-peak, while a symmetric 
two peak distribution appears at the X-point. Figure [7] 
shows the velocity distribution of ions with mj = 16m G 
at the same period. On the contrary to the electrons', a 
two peak distribution appears even at the O-point. The 
distribution implies that ions have not been thermalized 
yet due to the large mass. Figure [5] shows the electron 
velocity distribution at the same period for the case of 



3.6. The configuration of the magnetic field 

One can clearly see that the topology of the magnetic 
field changes when the filaments coalesce. In the frame- 
work of the resistive MHD, the coalescence of the self- 
confined plasma cy linders leads to the reconnection of 
the m agnetic field (|Parkerl [l957t iSweetl 119581 : iPetschekl 
1964), which is followed by the conversion of the mag- 
netic energy into th e kinetic ener gy of particles (see 
lUzdenskv fc Kulsrudll2000l : lKulsrudll2004n . However, re- 
sults of our simulation do not show such exchange of 
the energies. We have checked the velocity distribution 
around the X-point, where the magnetic reconnection 
might occur and found that it does not deviate from 
the Maxwell-Boltzmann distribution. This is at vari- 
ance with the prediction from the theory of the mag- 
netic reconnection. The lack of the release of the mag- 
netic energy seems to be attributed to the absence of 
the reconnection layer between the filaments, where the 
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Fig. 9. — T he c onfig urat ion of the stationary solution described 
by Equations (125 I t and H27I I. The color codes and arrows represent 
the density distribution of electrons (or ions) and the direction of 
the magnetic field in the x-y plane, respectively. 

anti-parallel magnetic fields coexist in a narrow and com- 
pressed region, does not form between the filaments. In- 
stead, magnetic fields vanish at the middle point between 
the filaments. 

4. QUASI-EQUILIBRIUM AND ITS TRANSITION 

From the analysis of results of the simulation, it is 
found that there is a quasi-equilibrium after the forma- 
tion of the filamentary structure in the saturation phase 
of the Weibel instability. In this section, we consider 
whether the equilibrium is expressed by an analytical so- 
lution of the Vlasov-Maxwell system. 

4.1. The analytical expression 

iSuzuki fc Shigevam a (2008) present a method to con- 
struct stationary solutions of the Vlasov-Maxwell system 
and derive a new two-dimensional equilibrium configu- 
ration of collisionless plasmas, whose velocity distribu- 
tion really resembles the quasi-equi librium shown above. 
Slight ly modifying the result of ISuzuki &: Shigevamal 
(2008)), we construct the following equilibrium, 



f e = C(v z - A z ) exp 



h = C(v z + A 2 ) 2 exp 




with the electromagnetic potentials expressed as 



4> = A X : 
which lead to 



0, A z = Bq sin a; sin y, 



(25) 
(26) 

(27) 



E x — E y — E z — 0, 



B x = B sin x cos y, B y = —B cos x sin y, B z = (118) 

One can easily check that these expressions satisfy the 
Vlasov-Maxwell system exactly. 

In the following we consider the case of mi = m e be- 
cause ions are not thermalized in the case of ui\ = 16m c . 
If the stationary solution described above actually corre- 
sponds to the quasi-equilibrium, free parameters in the 
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Fig. 10. — The v x , Hy-integrated electron velo city distr ibut ion of 
the stationary solution described by Equations (1 25 I t and (127ft . The 
solid and dashed lines represent the distributions at the O, and 
X-points, respectively. 

solution must be determined from results of the simula- 
tion. The value of t> (= Vi) is derived by fitting a gaussian 
to the x, y, v x , v z -integrated velocity distribution of elec- 
trons in the quasi-equilibrium phase. Then, we obtain 
v c = 0.14. The value of Bq is derived from the magnetic 
energy E m using Equations (f!2T)) and (j2"8|) as 

E m = -£- dx dy(sin 2 x cos 2 y + cos 2 x sin 2 £29) 

= n 2 Bl (30) 

On the other hand, E m ~ 0.2 from the result. So we 
obtain B ^ VO^/tt ~ 0.14. 

Figure [5] shows the density distribution and the config- 
uration of the magnetic field B x and B y of the equilib- 
rium described above. In Figure [TUI the solid and dashed 
line represent the v x , i^ -integrated velocity distributions 
given by Equations (|25|) and ([27]) at the O- and X-points 
with the parameters; v e = 0.14, and Bq = 0.14. The fila- 
mentary structure and the configuration of the magnetic 
field in Figure [5] remarkably resembles those in Figure 
[4j90 < t < 300). In addition, comparing Figure [TOl with 
Figure [3 one can see that the dominant features appear- 
ing in Figure [HJ the one-peak distribution at the O-point 
and the symmetric double-peak distribution at the X- 
point, are well reproduced by the analytical expressions 
([25]) and 1(27)1 . 

In the above discussion, we use the parameter v e and 
Bq determined from the result of our simulation. How- 
ever, the shape of the velocity distribution is sensitive to 
the parameters. The velocity distribution with different 
parameters, v e — 0.16 and Bq = 0.1, is shown in Figure 
llll which is more similar to the result of our simulation 
shown in Figure [8j 

Therefore, we conclude that the equilibrium configura- 
tion given by Equations |25)) and ((27)) provides an appro- 
priate expression of the quasi-equilibrium achieved in the 
beginning of the saturation phase, as long as the velocity 
distribution of electrons or that of ions with m\ = m e are 
concerned. 

4.2. Remaining problems 

Although we see that there is an analytical expression 
to describe the quasi-equilibrium configuration, some 
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problems remain. 

First problem is the disagreement of the distribution 
of ions with m\ — 16m e between the analytical expres- 
sion and the result of the simulation. As in Figure [3 
the velocity distribution of ions with m\ = 16m G have 
sharp peaks compared with the electron velocity dis- 
tribution, which cannot be reproduced by the expres- 
sions ([26)) and ([27)) . This disagreement arises since we 
truncate the series of v z exp(— v z /v?), (n — 0, 1,2,- • ■) 
at the second-order (n = 2) in constructing the veloc- 
ity distribution (|26|) . So this problem can be resolved 
by constructing higher-order stationary solutions of the 
Vl asov-Maxwell system accord ing to the procedure given 
bv lSuzuki fc Shigevamal (|2008f ); 

fj K [9j,o + 9j,xHx{v z /vj) + gj^Hzivz/vj) 

r 2 

+9j,3H 3 (v z /vj) H ] exp 



(31) 



where H n (v z /vj) are Hermite polynomials. Their coef- 
ficients gj t n are functions of A z . While this prescription 
enables us to make a more rigorous fits to the results of 
the simulation, the magnetic field cannot be expressed 
analytically. Furthermore, the expressions (|26[) and (|27[) 
are sufficient to reproduce the main feature of the result 
of the simulation, i.e., a single peak velocity distribu- 
tion at the O-point and a symmetric double peak veloc- 
ity distribution at the X-point. Therefore we take the 
construction of the higher-order stationary solutions as 
beyond the scope of the present work. 

Second is on the stability of the equilibrium. Results 
of the numerical simulation with the scale length of 2tt 
reveal that the quasi-equilibrium configuration achieved 
in the beginning of the saturation phase is unstable. The 
transition from the quasi-equilibrium to another appears 
in Figure [U This clearly indicates that the equilibrium 
configuration expressed by Equations ([25]) - ([28)) is unsta- 
ble. However, there is no analytical or numerical investi- 
gation into its stability and the long-term evolution. It is 
not clear from our simulation whether the critical scale 
length above which the equilibrium is stable exists or not. 
Nevertheless, we could make a brief speculation on the 
transition from the quasi-equilibrium state based on the 
result of the simulation. Comparing the time evolution of 
the kinetic energies K ex and K ez shown in Figure [31 and 
the time evolution of the filamentary structure shown in 
Figure [31 we can see a coincidence. At the beginning of 
the saturation phase (200 < t < 300), K ez takes greater 
values than K ex . Then, after the quasi-equilibrium con- 
figuration breaks, K ex and K ez oscillate changing their 
magnitude relation. Eventually, the amplitude of the 
oscillation decays and then K ex and K ez take similar 
values. This behavior coincides with the oscillation and 
the coalescence of the two filaments. In other words, 
the coalescence of filaments and the dissolution of the 
anisotropic distribution of the kinetic energies, K ex and 
K ezi occur simultaneously. 

5. DISCUSSIONS AND CONCLUSIONS 

In this paper, we have presented the behavior of a fila- 
mentary structure resulting from the Weibel instability. 
In order to simulate the behavior with high accuracy, we 
use a numerical scheme for the direct integration of the 
Vlasov-Maxwell system. The results reveal that there is 
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Fig. 11.— Same as Figure ITOl but for v e = 0.16 amd B = 0.1. 



a quasi-equilibrium configuration after the saturation of 
the Weibel instability. Then, it turns into another config- 
uration. The quasi-equilibrium configuration is found to 
be described by the analytical expression ([25 |) - ([28|) . The 
coalescence of the filaments is synchronized with the time 
evolution of the kinetic energies of electrons K ex and K ez . 

Comparing results of the simulation, we con- 
sider the validity of mo dels proposed so far. 
iMilosavlievic fc Naked (|2006f ) treated the filamen- 
tary structure in the framework of the MHD. In our 
simulation, however, the velocity distribution at the X- 
point is far from the Maxwell-Boltzmann distribution as 
seen in Figures [51 and [51 which implies the treatment 
of the plasma as a magnetized fluid is not appropriate. 
Furthermore, although they assumed that the filament 
is isolated and other filaments assert a negligible effect, 
two peaks in the velocity distribution at the X-point 
resultant from our simulation show that the outskirts 
of filaments coexist around t he X- point. The model 
proposed by iMedvedev et al.1 (|2005l ) is based on a toy 
model. Although this model ignores the effect of the 
dissipation of the magnetic field energy via reconnection 
of the magnetic fields, it is expected in the framework of 
the resistive MHD. Our simulation docs not show any 
sign of the reconnection. As stated earlier, we attribute 
the lack of the magnetic reconnection to the absence 
of the reconnection layer between the filaments. But, 
this can also be due to the coarseness of the simula- 
tion. While judging whether the lack of the magnetic 
reconnection is real or not needs careful investigation, 
the results of our simulation just i fy the treatment of 
the co alescence in IMedvedev et al.l ()2005l ) . IChang et al.l 
(2008) modeled the decay of the magnetic turbulence 
generated by the Weibel instability. This treatment 
is appropriate only in the region where the magnetic 
trappi ng is not so important. Actually, Chan get al.l 
(2008) confirm that their model cannot predict the time 
evolution of magnetic field strength in the presence of a 
magnetic field with long wavelength. 

While the detailed dynamical behavior of the filamen- 
tary structure is revealed in this study, some problems, 
for example, the long-term evolution of the filaments and 
whether the critical scale-length above which the equilib- 
rium is stable exists or not, remain unclear. In addition, 
because the distribution function is assumed to be uni- 
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form along the direction of the initial bulk flow of plas- 
mas, we cannot deal with some 3D perturbation on the 
equilibrium, for example, bending and twisting of the fil- 
aments. The stability of the equilibrium described by 
Equations (125 |) - (|28|) . which is presumably a key ingredi- 
ent for a strong magnetic field to survive in long term 
and serve as an accelerator of particles in the transition 
layer of a collisionless shock front, should be investigated 
rigorously including 3D effects. 
Although this study treats non-relativistic plasmas, 



relativistic effects become vital in some astrophysi- 
cal phe n omen a, e.g., GRB afterglows, jets from AGN. 
I Suzukil (l2008f) prov i des t he relativistic extension of 
Suzuki & Shig evamal (|2U 1)M. i.e, stationary solutions of 
the relativistic Vlasov-Maxwell system. It is not irrel- 
evant that on e expects the equilibrium constructed by 
ISuzukil (|2008f) to provide an appropriate model for the fil- 
amentary structure resulting from the Weibel instability 
in relativistic plasmas as the non-relativistic counterpart 
does. 
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